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ABSTRACT 



Among the hot Jupitcrs that transit their parent stars known to date, the two 
best candidates to be observed with transmission spectroscopy in the mid-infrared 
(MIR) are HD 189733b and HD209458b, due to their combined characteristics of 
planetary density, orbital parameters and parent star distance and brightness. 
Here we simulate transmission spectra of these two planets during their primary 
eclipse in the MIR, and we present sensitivity studies of the spectra to the changes 
of atmospheric thermal properties, molecular abundances and C/0 ratios. Our 
model predicts that the dominant species absorbing in the MIR on hot Jupiters 
are water vapor and carbon monoxide, and their relative abundances are deter- 
mined by the C/0 ratio. Since the temperature profile plays a secondary role 
in the transmission spectra of hot Jupiters compared to molecular abundances, 
future primary eclipse observations in the MIR of those objects might give an 
insight on EGP atmospheric chemistry. We find here that the absorption features 
caused by water vapor and carbon monoxide in a cloud-free atmosphere, are deep 
enough to be observable by the present and future generation of space-based ob- 
servatories, such as Spitzer Space Telescope and James Webb Space Telescope. 
We discuss our results in light of the capabihties of these telescopes. 

Subject headings: Atmospheric Effects, Occultations, Radiative Transfer, Techniques: 
spectroscopic 
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Introduction 



Extrasolar Giant Planets (EG Ps) are now being discovered at an accelerating pace 
( jSchneiderl l2006l : iButler et al.ll2006l ). In particular, an increasing interest has been focused 
on hot Jupiters that transit their parent stars, since they represent a valuable tool to 
determine key physical parameters of the EGPs, such as atmospheric composition a nd 
dynamics, thermal properties and presence of condensates (jSeager and Sasselovl |2000| ) . 



The most studied transiting extrasolar planet, HD209458b, orbits a main sequence 
G-type star at 0.046 AU (period 3.52 days). It is the first one f or which rep e ated 
transits across the s tellar disk were observed (~ 1.6% absorption; iHenrv et al.l (120001): 
Charbonneau et al. ( 2000l )). Along with radial velocity measurements ( Mazeh et al. 2000l ). 



it was possible to determine mass and radius (Mp ~ 0.69 Mj^p, Rp ~ lARjup), confirming 
the planet is a gas giant, with one of the lowest densities discovered so far. Owing to this 
property, its very extended atmosphere is one of the best candidates to be probed with 
transit techniques. In particular the upper atmosphere extends beyond the Roche lobe, 
showing a population of escaping atoms. This discovery was possibly due to the observed 
extraordinary deep absorpti ons in HI, 01 and CII ove r the stellar emissions lines (15%, 
13% and 7.5% respectively; IVidal-Madjar et al.l (120031 . 120041 )). The numerous follow-up 
observations of HD209458b, also in clude the detection and upper limits of absorptio n 
features in the deeper atmosphere (Charbonneau et al. 20021: [Richardson et al. 2003a| jbl: 
Deming et al.ll2005bl : iRichardson et al.ll2006l ). Most recently iDeming et al.l (l2005al ) detected 
the thermal emission of this p lanet with the Spi t zer S pace Telescope during a secondary 
transit in the 24 fim band and IRichardson et al.l (120061 ) detected the first primary eclipse in 
the same band. 

To explain these observations several models were proposed, including atmospheric 
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so far. Additional observations are required, though, to constrain the past, present and 
future modeling effort. 



The planet HD189733b, recently discovered by iBouchy et al.l (|2005l ). with mass 
^ 1.15 Mj„p and Rp ~ 1.26 i^j^p, orbits an early main sequence K star at 0.0313 A 
is an exoplanet transiting the brightest and closest star discovered so far. 



Mp ~ 



Here we focus our interest on both planets HD209458b and HD189733b, and we 
model the spectral absorption features in the Mid-Infrared (MIR) due to the most 
abundant atmospheric molecules during their primary eclipse, i.e., when the planet passes 
in front of the parent star. The use of transmission spectroscopy to probe the upper 
layers of the tra nsiting EGPs, has been particularly successful in t he UV and visible 
spectral ranges ( Charbonneau et al. 2002; Richardson et al. 2003a| jbl: IPeming et al. 2005b; 
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Vidal-Madjar et al.l |2003| . |200J) , and only very recently attem pted in the MIR, in the 24 



/im band, using Spitzer observations (IRichardson et al.l 120061 ). Two circumstances make 
it an approach worth considering. First, the high surface temperatures, relatively small 
masses, and mostly hydrogen atmospheres of close-in EGPs imply large atmospheric scale 
heights. As a consequence, for spectral features that span a reasonably wide bandwidth, so 
that the total photon flux is not too small, this is a feasible and diagnostically powerful 
technique. Second, this is a complementary approach to the secondary eclipse observations. 
Transmission spectroscopy, taken during primary eclipse, is sensitive to different parameters 
and regions of the atmosphere, compared to emission spectroscopy, on which the secondary 
eclipse method is based. 



Oxygen versus Carbon H2O versus CO In a solar system like ours, a significant 
amount of water vapor (H2O) can exist only in planetary atmospheres at orbital distances 
less than 1 AU. The requirement is certainly met for the known transiting EGPs. Carbon 
monoxide (CO) and methane (CH4), and other photochemical products, such as carbon 
dioxide (CO2) and acetylene (C2H2), are plausibly present in the atmospheres of EGPs, 
and possibly abundant to be detected. These species have strong absorption bands in 
the MIR, and more importantly, in spectral regions compatible with present and future 
space-base d observations such as the Spitzer Space Telescope or the James Webb Space 
Telescope (IGardner et al.ll2006l ). Given O and C, H2O and CO will be controlled mainly by 
the relative abundances of these two species. 



If C/0 ratio is close to the solar, H2O, CO and CII4 abundances a re determined by 
the t hermodynamic equilibrium chemistry in the deep atmosphere (ILiang et al.ll2003l . 



2004f ). 



If C/0 ratio is above solar, the atmospheric chemistry might change d ramatically 
and, according to the scenario proposed by iKuchner and Seagerl (120061 ) . planets 
should show a significant paucity of water vapor in their atmospheres, carbon rich 
species should by contrast be enhanced. In particular, CO is expected to be the 
dominant carbon- bearing molecule at high t emperatures and CH4 the dominant at 
low temperatures (IKuchner and Seagerl |2006| ). 



• If C/0 ratio is below solar, the atmosphere is depopulated of carbon-bearing molecules 
and water-vapor is the dominant species between ~ 10~^° to 1.5 bars. 



2. Description of the Model 



We have built a model of planetary atmosphere and calculated the expected absorption 
of the stellar light when filtered through the planetary atmospheric layers. This has 
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been already discussed in the literature. In particular for o ur simulations, we have used 
the geometry and the e quations described in Brown f 200ll ) (fig. 1, configuration 2) and 
Ehrenreich et al.l (120061 ) (sec. 2.1, fig. 1). Our cloud and haze-free atmospheres were divided 



in forty layers spanning from ~ 10 to 1 bars. 

Photochemical models are used to determine the molecular abundances of 33 species 
above ~1 bar altitude level. We start with four parent molecules H2, CO, H2O, and CH4. 
Their relative abundances are determined by thermochemistry in the deep atmosphere, and 
are fixed as our lower b o undary coii dition. Chemical reactions and eddy mixii ig profiles ar e 
taken from 



Liang et al.l (l2003l . |200J). Details of the model can be referred to iLiang et al 



( 120031 . 12OO4I ) and references contained therein. For the simulation of the photochemistry 
of HD209458b, we adopt the solar spectrum. For HD189733b we use the spectrum of 
HD22049, which is a K2V star similar to HD189733 (ISegura et al.ll2003[ ). We have repeated 
our calculations for three temperature-pressure profiles (fig. [T] right), to test the sensitivity 
of our results to these assumptio ns. The modeled chemical abundances show a negligible 



dependence on temperature (see (iLiang et al.ll2004 )) 



The a bsorption coefficients in th e MIR were estimated using a line-by-line model, 
LBLABC ( iMeadows and Crispl Il996l ). that gen erates monochr o matic gas abso rption 



coefficients from molecular line lists -HITEMP, iRothman et al.l (lin preparation 
of the gases present in the atmosphere. 



for each 



3. Results 

Fig. [1] (left) shows the molecular profiles of H2O, CO, CH4, CO2 and C2H2 for both 
planets HD189733b and HD209458b, calculated by the photochemistry model with solar 
C/0 ratio as boundary condition. On HD189733b, H2O, CH4 and C2H2 are more abundant 
in the upper atmosphere compared to HD209458b, since HD189733 is a later type star, 
therefore the photo-dissociation processes occurring in the atmosphere of that planet are 
less significant. 

Sensitivity to molecular abundances and C/0 ratio. Figs. [2] show the predicted 
absorption signatures due to water vapor and CO on the planets HD189733b and 
HD209458b. The three plots compare the spectral absorptions of these two species when 
C/0 ratio is solar (standard case, solid line, see fig. [1] left for mixing ratios) and when 
is below and above solar. As specific examples, we have assumed H2O to be 10 times 
more and at the same time CO 10 times less abundant than the standard case (dotted 
line), and vice- versa (dashed line). When we increase/diminish H2O and CO of a factor 
10, the absorption is increased/decreased by a constant ~ 0.03 % through all the selected 
wavelength range. 



In these figures, we show also the signatures (white rhombi, squares and triangles) 
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relative to the three cases (standard, C/0 ratio below and above solar) averaged over 
the IRAC, IRS and MIPS bandpasses (centered at 3.6, 4.5, 5.8, 8, 16 and 24 fxra), 
the instruments on board the Spitzer Space Telescope (see table 1 for the calculated 
absorptions). Water vapor has strong absorption lines through all the selected spectral 
range in the MIR, the CO signature appears only in a narrower wavelength interval, where 
the IRAC 4.5 fim bandpass is centered. When the C/0 ratio is above solar, the triangle is 
expected to appear above the rhombus in that band, indicating the strong CO contribution 
to the total absorption (table 1, numbers in bold). 

Note that we included the species CH4, CO2 and C2H2 in our calculations, which 
strongly absorb in the MIR. However their abundances are too small compared to CO 
and H2O, and they are masked by these species. For example, we show in fig. [3] (left) 
the contribution due to CO2 in the case of C/0 above solar ratio. CO2 is inc reased by a 



facto r 10 here, to follow CO behavior consistently with chemistry predictions (iLiang et al.l 



20031 ) . When CO2 is present, an increase in absorption of ~ 0.02 % is found in the band 
centered on 15 /Um and a very narrow peak reaching 0.07 % is visible at shorter wavelengths. 
Although not negligible, the contribution due to CO2 is masked by CO and water. 

Sensitivity to temperature. The effects due to temperature variations (fig. [3] right) 
are negligible compared to the changes in molecular mixing ratios. Temperature plays a 
secondary role in the determination of the optic al depth: it affects the ab sorption coefficients 



and the atmospheric scale heights (see eq. 2 in lEhrenreich et al.l (120061 )). For HD189733b, 
the discrepancy between the standard and the hot profile is less than 0.005%, and between 
the standard and the very hot profile has a maximum of 0.014% in the 15-30 //m range (fig. 
in] right). Analogous results are obtained for HD209458b. 



Discussion 



I n our model we did not include the contri bution of hazes or clouds (lAckerman and Marley 



2OOII : iLunine et al.lll989l : iFortney et al.ll2005l ). Due to their presence, the atmospheric 



optical depth might increase, partially masking the absorption features due to atmospheric 
molecules. In the case of water vapor and CO, only clouds/hazes lying at altitudes higher 
than 1 bar might affect our results. Predictions of cloud/haze is particularly difficult for 
EGPs, since the few observations we have are not sufficient to constrain all the cloud 
microphysics and aerosol parameters. Moreover, the planetary limb observable during 
the star occultation might show the signatur es of both the night and day side of those 
planets, which are presumably tidally locked (jiro et al.l 120051 ). The thermal profiles, hence 
the condensate dynamics, might be very different on the two sides. Consequently, a 
more complete model able to predict cloud and haze location and optical characteristics, 
should contain a 3-D dynamical simulation of the atmosphere. In this paper we limit our 
simulations to the cloud-haze free atmosphere, with the caveat they might be perturbed by 
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the possible presence (constant or variable) of optically thick particles in the atmosphere 
above 1 bar pressure. 

Same considerations are valid for the thermal profiles. A extensive literature is available 
on T — P profiles for EGPs at pressures from ~1 bar to 10~^-10~^ bars, most recently 



including 3D dynam ical effects (jShowman and Guillotll2002l : liro et al.ll2005l : ICho et al.ll2006 



Burrows et al.l 120061 ). For transmission spectroscopy in the MIR, we need to consider als o 
the co n tribut i on of the upper atmosphere. The T — P profiles calculated by iTian et al. 



(I2OO5I ): lYelld (I2OO4I ) suggest that the trend for the atmospheric temperature is to increase 
in the exosphere. For our simulations, we use a. T — P profile compatible with the lower 
atmosphere models cited above up to 10~'^-10'^ bars, and then we consider three cases: 
the atmospheric temperature decreases up to 10~^° bars (standard), the atmosphere is 
isothermal (hot, very hot profiles). Our results, show that the differences among the spectra 
calculated with the three profiles are within 0.009 % for A < 14/im and within 0.014 % at 
longer wavelengths, so we are confident our simulations will not significantly change using a 
more refined thermal structure. 

Our model atmospheres extend to 10 ~^° bars, where non local thermodynamic 
equilibrium (non-LTE) effects might occur (IKutepov et al.l 119981 ). However, if we truncate 
our calculations to 10~^ bars, we obtain a maximum error of ~ 0.02% at 30 fim (no 
discrepancy for wavelengths shorter than 20 ^m), indicating that our calculated absorptions 
in LTE regime are correct at first order approximation. 

In order to detect the presence of water vapor and CO on HD 189733b and HD209458b 
in the MIR, an extra absorption of ~ 0.15 % is expected to be added to the 2.85 % 
and 1.6 % due to the optically thick disks at 1 bar atmospheric level. To estimate the 
chemical abundances, an accuracy of at least 0.03 % is needed. By inspection of the relative 
absorption of the IRAC 4.5 /im bandpass with respect to the others, we might be able to 
infer the C/0 ratio, but in this case an extremely high S/N is required. HD189733 is a 
bright KOV star of magnitude K = 5.5. We estimate the brightness in the four IRAC bands 
to be of the order of 1850, 1100, 730 and 400 mJy at 3.6, 4.5, 5.8 and 8 /im respectively. 
For HD209458, a GOV star with K-magnitude of 6.3, the IRAC predicted fluxes are 878, 
556, 351 and 18 9 mJy. According t o these numbers, a better S/N should be obtainable 
for HD189733b (jPeming et al.ll2006l ). and this makes HD189733b a better candidate for 
observations. 



5. Conclusions 

In this paper we have presented simulations of transmission spectra of two extrasolar 
giant planets during their transit in front of their parent star. According to our calculations, 
we estimate an excess absorption in the IR of up to 0.15 % for HD189733b and up to 
0.12 % for HD209458b (C/0 ratio - solar), in addition to the nominal 2.85 % and 1.6 % 
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absorptions measured at shorter wavelengths. If water were far less abundant, other species 
might be observable, depending on their mixing ratios. Among them, CO2, CH4 and C2H2 
are the best candidates. 

According to our simulations, transmission spectra of EGPs in the MIR arc sensitive 
to molecular abundances and less to temperature variations. Temperature influences the 
transmission spectrum by way of its influence on the atmospheric scale height, as discussed 
by Brown (2001), and on the absorption coefficients. 

If water vapor and CO are as abundant as photochemical models predict, we expect 
they can be detected with the IRAC, IRS and MIPS instruments on board the Spitzer 
Telescope and with future telescopes hke JWST. Moreover, if an accuracy of 0.03 % is 
obtainable, future observations may give a flrst direct estimate of H2O and CO abundances 
in the upper atmosphere of EGPs and possibly -depending on their mixing ratios- a 
constraint on CO2, CH4 and C2H2. 



Acknowledgments 

We would like to thank the anonymous referee for his help to improve the paper, L. S. 
Rothman for having provided the HITEMP data list, R. Ferlet, J. M. Desert, F. Bouchy, G. 
Hebrard, A. Noriega Crespo and S. Carey, for their valuable inputs, and C. D. Parkinson 
for useful comments. G. Tinetti is currently sponsored by the European Space Agency. M. 
C. Liang and Y. L. Yung are supported by the NASA grant NASA5-13296 to the California 
Institute of Technology. 



-9- 



REFERENCES 

Ackerman A.S. and Marley M. S. (2001), ApJ, 556, 872. 

Bouchy F., Udry S., Mayor M., Moutou C, et al. (2005), A&A, 444, L15. 

Brown T. M. (2001), ApJ, 553, p 1006. 

Burrows A., Sudarsky, D., and Hubbard, W. B., (2003), ApJ., 594, 545. 
Burrows A., Sudarsky D. and Hubeny I. (2006), ApJ., accepted, astroph 0607014. 
Butler R. P., J. T. Wright, G. W. Marcy, et al. (2006), ApJ, 646, 505. 
Charbonneau D., Brown T. M., Latham D. W. and Mayor M. (2000), ApJ, 529, L45. 
Charbonneau D., Brown T. M., Noyes R. W. and GiUiland R. L. (2002), ApJ 568, 377. 
Cho J. Y. K., Menou K., Hansen B. M. S. and Seager S. (2006), astroph 0607338. 
Deming D., Seager S., Richardson L. J. and Harrington J., (2005a), Nature, 434, 740. 
Deming D., Brown T. M., Charbonneau D. et al. (2005b), ApJ, 622, 1149. 
Deming D., Harrington J., Seager S. and Richardson L. J., ApJ, 644, 560-564 
Ehrenreich D., Tinetti G., Lecavelier A. et al. (2006), A&A, 448, 379. 
Fortney J. J., Sudarsky D., Hubeny 1, et al. (2003), ApJ, 589, 615. 

Fortney J. J., Marley M. S., Lodders K., Saumon D. and Freedman R. (2005), ApJ, 627, 
L69. 

Gardner J. P., Mather J. C.; Clampin M. et al.. Space Science Reviews, in press. 

Henry G. W., Marcy G. W., Butler R. P. and Vogt S.S. (2000), ApJ, 529, L41. 

Iro N., Bezard B., Guillot T. (2005), A&A, 436, 719. 

Kuchner M. and Seager S. (2006), astroph 0504214, ApJ, submitted. 

Kutepov A. A., Gusev O. A., and Ogibalov V. P., 1998, J. Quant. Spectrosc. Radiativ. 
Transfer, 60, 199-220. 

Lecaveher des Etangs A., Vidal-Madjar A., McConnell J. C. and Hebrard G., et al. (2004), 
A&A, 418, LI. 

Liang M. C., Parkinson C. D., Lee A. Y.-T. et al. (2003), ApJ, 596, L247. 



- 10 - 



Liang M. C, Seager S., Parkinson C. D., et al (2004), ApJ, 605, L61. 

Lunine J. L, Hubbard W. B., Burrows A., et al. (1989), ApJ, 338, 314. 

Mazeh T., Naef D., Torres G., et al. (2000), ApJ, 532, L55. 

Meadows V. S. and Crisp D. (1996), JGR, 101, 4595. 

Richardson L. J., Deming D. and Seager S., (2003a), ApJ, 597, 581-589. 

Richardson L. J., Deming D., Wiedemann G., et al. (2003b), ApJ, 584, 1053. 

Richardson L. J., Harrington J., Seager S. and Deming D., (2006), ApJ, 649, 1043-1047. 

Rothman L. S., C. Camy-Peyret, J.-M. Flaud, et al., HITEMP, the High- Temperature 

Molecular Spectroscopic Database is being prepared for the Journal of Quantitative 
Spectroscopy and Radiative Transfer. 



Schneider J. (2006), The Extrasolar Planets Encyclopaedia, http://exoplanet.eu/ 
Seager S. and Sasselov D. D. (2000), ApJ, 537, 916. 

Seager S., L. J. Richardson , B. M. S. Hansen et al. (2005), ApJ, 632, 1122. 

Segura A., Krelove K., Kasting J. P., et al. (2003), Astrobiology, 3(4), 689 -708 

Showman, A. P. and Guillot T., (2002), A&A, 385, 166. 

Tian P., Toon O. B., Pavlov A. A. and De Sterck H. (2005), ApJ, 621, 1060. 

Vidal-Madjar A., Lecavelier des Etangs A., Desert J. M., et al. (2003), Nature 422, 143. 

Vidal-Madjar A., Desert J.-M.; Lecavelier des Etangs A. et al.(2004), ApJ, 604, L69. 

Yelle R. (2004), Icarus 170, 167. 



This manuscript was prepared with the AAS L^T[t;X macros v5.2. 



- 11 - 



Table 1. Calculated absorptions averaged over IRAC, IRS & MIPS bandpasses.We recall 
that 2.85 % and 1.6 % are the nominal absorptions due to the optically thick disks at 1 bar 
atmospheric level. In bold the absorptions when C/0 ratio is above solar and CO strongly 

contributes in the 4.5 /xm IRAC band. 





C/0 ratio 


3.6 fxm 


4.5 /im 


5.8 /im 


8 /im 


16 /im 


24 /im 


HD189733b 


< solar 


2.953 


2.970 


3.004 


3.001 


3.006 


3.024 


HD189733b 


solar 


2.930 


2.959 


2.978 


2.974 


2.977 


2.992 


HD189733b 


> solar 


2.905 


2.961 


2.960 


2.952 


2.957 


2.973 


HD209458b 


< solar 


1.715 


1.734 


1.773 


1.770 


1.776 


1.796 


HD209458b 


solar 


1.690 


1.722 


1.746 


1.741 


1.747 


1.766 


HD209458b 


> solar 


1.663 


1.725 


1.724 


1.714 


1.720 


1.738 




Fig. 1.— Left: Profiles of H2O (solid line), CO (dotted line), CO2 (dot-dashed line), CH4 (dashed line), 
and C2H2 (triple-dot-dashed line) for planets HD209458b (black lines) and HD189733b calculated with the 
photochemistry model described in iLiang et al. I (|2003ll2004 . Right: temperature-pressure profiles used for 
our simulations. Solid line: standard, dotted line: hot profile, dashed line: very hot profile. 
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Fig. 2. — Modeled spectral absorptions of water vapor and CO in the atmospheres of HD189733b (left) and 
HD209458b (right) during their transits. Solid line: standard case, C/0 ratio — solar (fig. [l]left); dashed 
line: CO is over- and H2O under-abundant of a factor 10; dotted line: the opposite. White triangles, squares 
and rhombi indicate the spectral absorptions averaged over the Spitzer's IRAC, IRS and MIPS bandpasses, 
centered at 3.6, 4.5, 5.8, 8, 16 and 24 //m. Only in the 4.5 /im band, where CO shows a strong absorption, 
the triangle and rhombus overlap. 




Fig. 3. — Left: residue spectrum from difference of 2 curves: dashed line plot in fig. [2] top with and 
without the contribution of CO2. CO2 is assumed to be, like CO, over-abundant by a factor 10. Right: 
residue spectra from difference of 2 curves. Dotted line: difference from hot and standard profiles; dashed 
line: difference from very hot and standard profiles (fig [1] right). The discrepancies amidst the curves is 
negligible compared to the changes in molecular abundances (fig. [2]). 



